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ABSTRACT 


Theoretical rate constants are presented for excitation 
of NH 3 by collisions with He. The lowest 22 levels of ortho-NH 
and the lowest 16 levels of para-NH 3 are considered at kinetic 
temperatures of 15 -■» 300 K. 
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I. Introduction 

Analysis of microwave and IP. observations of interstellar 
ammonia requires knowledge of the rate constants for collisional 
excitation among the observed rotation-inversion levels. It 
is not currently possible to obtain the required state-to- 
state rate constants for excitation by H 2 and He experimentally. 
However, some related information is available; in particular, 
pressure broadening data provide a measure of the total rate 
of excitation out of the two spectral levels, and microwave 
double resonance data provide a measure of the relative 
rates among different levels. On the other hand, these 
collision rates can be obtained theoretically by solving the 
electronic Schrodinger equation for the intermolecular forces, 
followed by molecular scattering calculations to obtain 
collision cross sections. Such calculations can also predict 
pressure broadening and double resonance results so that 
comparisons with available experimental data provide some 
check on the theory. 

A theoretical study of NH^ excitation by collisions with 
He was recently presented. ^ In that study the intermolecular 
forces were obtained by combining a previous Hartree-Fock 
calculation for the short-range forces with estimates for the 
long-range forces. Scattering calculations were performed 
for para-NHg-He to obtain state-to-state cross sections and 
rate constants and also to predict pressure broadening and 
microwave double resonance results. Agreement with pressure 


3 


broadening data was excellent. Agreement with double 
resonance data, however, was mixed. It was argued that the 
discrepancies here are due to an oversimplified analysis of 
the experimental data and do not indicate a serious error 
in the theoretical calculations. 

Due to molecular symmetry and nuclear spin, ammonia 
molecules are either para-NH 3 for rotation-inversion levels 
with K-0 , 3 , . . . , 3n or ortho-NHg for levels with K-l, 2, 4, 5 , . , . , 
2n+l. These are interconverted only very slowly by radiation 
or by non-reactive collisions and, for all practical purposes, 
can be considered as distinct species under the conditions 
in interstellar space. Ref. 1 presented calculations only 
for para-NH^. The present work extends these studies to 
obtain collision rates for ortho-NHg. For convenience, rates 
for para-NH 3 are also tabulated here for more kinetic temper- 
atures than given in Ref. 1. Because the computational 
methods are essentially the same as those used previously, 
only a brief description will be given in the following Section. 
The accuracy of the resulting rate constants is addressed 
in Section III, 

II. Computational Details 

The interaction labeled ’"potential 1" in Ref, 1 was used 
here. The short-range forces are from Rartree-Fock calcula- 
tions for NH 3 -He. To this are added estimates for the 
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following long-range forces: a spherical* R dispersion term, 

the quadrupole-induced-dipole interaction, and the anisotropy 

in the dispersion. As noted in Ref. 1, the value used for 

this last term is probably too large by about a factor of four, 

but it has only a minor effect on numerical results. 

Scattering calculations were done within the coupled 

states (CS) approximation. The accuracy of this approximation 

has been tested by comparisons with essentially exact close 

1 2 

coupling calculations. ' In general, CS results for the 
larger, more important cross sections are accurate to better 
than 10%; even for less probable transitions the error is 
rarely more than 50%. The NH 3 rotational wave functions are 
taken as linear combinations of symmetric top functions with 
quantum numbers J and K: 

| JK e > - I JK> + € I J-K>. 

. *T 

For K=0, only e=+l is permitted; for K>0, e-+(-l) where the 

upper sign refers to the lower member of an inversion doublet 

and vice versa . Rotational energies were computed from the 

standard formula for symmetric tops with rotation constants 
-1 -1 

A=B=9.9402 cm and 0=6.3044 cm . The separation between 
members of an inversion doublet (about 0.8 cm” 1 ) was ignored 
in calculations for ortho-NH^. The effect of these splittings 
was considered in the earlier calculations for para-NH 3 and 
found to be unimportant. 
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Scattering calculations were done at several collision 

energies. Because the astrophysical sources generally have 

low kinetic temperatures, special attention was given to 

energies near excitation thresholds. The collisional reduced 

mass was taken as 3.24 atomic mass units. Enough energetically 

closed channels and enough partial waves were included to 

assure convergence of inelastic cross sections to better than 

a few percent. In particular, calculations were done including 

the 12 rotation-inversion levels with J*0-5 for nine energies 

-1 

from 25 to 145 cm ; including the 17 levels with J*Q~6 for 

-1 

ten energies from 160 to 325 cm ; including the 22 levels 

—1 

with J=0-7 for six energies from 360 to t>00 cm ; and including 

the 27 levels with J«0~8 for energies of 800, 1000, and 1200 
-1 

cm , Resulting cross sections were integrated over Boltzmann 
distributions of collision velocity to obtain the rates 
presented in Table 1, 

The cross sections obtained previously for para-NH^ 
have been integrated over Boltzmann distributions of collision 
energies to obtain the rates given in Table 2, 

III. Accuracy of Results 

As discussed in Ref. 1, pressure broadening provides 
a measure of the total inelastic collision rate, and the 
calculations presented there agreed with microwave data to 
about 10%. In addition, broadening of the pure rotational 
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-1 

J,K * 0,0 -1,0 infrared line at about 20 cm has been 

measured 3 giving an effective nross section of 17 R 2 . The 
present calculation tor vhis quantity of 18 R 2 is in good 
agreement. 

Agreement of theory with microwave double resonance 
data for para-NH 3 was found to be mixed. Fewer data are 
available for ortho-NH 3 ; however, comparison with the 
present theoretical results shows the same qualitative 
features found for para-NH 3> As discussed in Ref. 1 this 
lack of accord is believed to reflect an oversimplification 
in the analysis of the data and not an inherent flaw in 
the theoretical collision rates. 

The major source of error in the present results is 
undoubtedly due to uncertainties in the intermolecular 
forces. Ref. 1 considered the effect on collision cross sections 
of reasonable variations in the long-range interaction.. 

It was found that the larger cross sections were relatively 
insensitive to variations in the potential, experiencing 
changes of a factor of two or less, suggesting that the 
present results are at least this accurate. 

Thus, from comparison^ with pressure broadening data, 
the total excitation rate Out of a given level is likely to 
be predicted within 20% of the correct value in the present 
study. The absolute values of individual state- to- state 
rates are probably accurate to better than 50%, although 
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somewhat larger errors, especially in the smaller, less 
important rates, cannot be ruled out with complete confidence. 
Also, as discussed in Kef. 1, values for para-NH 3 in Table 2 
required extensive extrapolation in the high energy tail of 
the Boltzmann distribution for temperatures above about 
200 K; however, this is not expected to introduce errors 
larger than those due to uncertainty in the potential. 

Excitation by collisions with cold H 2 (in its lowest, 

J-0 rotational level) is expected to be similar to excitation 
by He, and this is supported by pressure broadening data. 1 
Owing to lower mass, at a given temperature H 2 molecules 
have a 40% higher velocity than He atoms. To account for 


-1~ V* 4 M 4k 44 -3 r n 44. JL» 44. — • *1 m jM U. X- *1 . . 1 44 M 44 4k l— 4a. -4. 4-4 4 4 4k -J W M. 4k. «. * 4 44 « Lk 44k 4k 4 J 4 U 4 »4 4> 

uu,io ciiua j.ui ums: aia.yui-xj xaitjcx uwocivcu pj. coo uj. c uxuaucnxuy 


cross sections, it is recommended that the rates in Tables 
1 and 2 be increased by 50% for excitation by cold H 2 . If 
the H 2 molecules are in higher rotational levels (J>0) they 
are more effective in exciting NH^ due to the long-range 
quadrupole-dipole interaction. This enhances the rate of 
dipole allowed transitions in NH 3 , i.e., AJ^Q transitions 
across an inversion doublet and AJ-1 transitions from lower- 
to-lower and upper-to-upper members of inversion doublets. 
Although detailed calculations are not available for this 
system, enhancement is probably about a factor of two for 
these cases. This is in harmony with the observation that 
the total excitation rate, as measured by pressure broadening, 
is about 50% greater for H 2 (J>0) than for H 2 (J-0). 
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